Low temperature spin diffusion in the 
one-dimensional quantum 0(3) nonlinear cr-model 
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An effective, low temperature, classical model for spin transport in the one-dimensional, gapped, 
quantum 0(3) non-linear cr-model is developed. Its correlators are obtained by a mapping to a 
model solved earlier by Jepsen. We obtain universal functions for the ballistic-to-diffusive crossover 
and the value of the spin diffusion constant, and these are claimed to be exact at low temperatures. 
Implications for experiments on one-dimensional insulators with a spin gap are noted. 



Over the past decade, a large number of one- 
dimensional, insulating, Heisenberg antiferromagnets 
with a zero temperature (T) spin gan-have been stud- 
ied: these include integer spin chainsoa and half-integer 
spin-ladder systemsu. In the large spin S limit, thcJow 
energy properties of these compounds are describedu by 
the one-dimensional quantum 0(3) non-linear cr-model 
(without any topological term), and there is evidenceu 
that the mapping to this continuum model is quantita- 
tively accurate even for the S — 1 spin chain. Theo- 
retically, muchisinown about the quantum field theory 
of the cr-modelETH, and this information has been valu- 
able in understanding the properties of the spin chains. 
The low energy spectrum of the cr-model consists of a 
triplet of massive particles, and their ballistic propaga- 
tion describes many exactly known dynamic correlations 
at T = 0. For T > however, exact results have so far 
been limited to static, thermodynamic observablesa. 

In this paper, we obtain dynamic, non-zero T correla- 
tors using a semiclassical methodcH: we claim that all of 
our results are asymptotically exact at low T, but this 
has not been rigorously established. We present univer- 
sal functions which describe the crossover from ballistic 
spin transport at short scales, to diffusive behavior at 
the longest scales; as a bi-product, these functions yield 
the exact value of the spin diffusion constant. The na- 
ture of spin transport for any small T > is therefore 
qualitatively different from that at T = 0. 

The imaginary time (r) action of the cr-model is 



c f 1/T 
A = — / drdx 



(d x n a ) 2 + —(d T n a - ie aj3 ~,H pn^f 



where x is the spatial co-ordinate, a,/3, 7 = 1,2,3 are 
0(3) vector indices over which there is an implied sum- 
mation, e Q/ 3 7 is the totally antisymmetric tensor, c is a 
velocity, H a is an external magnetic field, and the par- 
tition function is obtained by integrating over the unit 
vector vield n Q (x,r), with tv^{x,t) = 1. We use units 
in which fi = ks = 1 and have absorbed a factor of the 
electronic magnetic moment, g e fJ-B, into the definition 
of the field H. The dimensionless coupling constant g 
is determined by the underlying lattice antiferromagnet 



at the momentum scale A ~ inverse lattice spacing to 
be g ~ 1/5. We shall only be interested in the physics 
at length scales ^> A -1 and time scales 3> (cA) , this 
physics is universally characterized by the dimcnsionful 
parameters c, H, T, and A, the energy gap at T = H = 0. 
The magnitude of A is determined by non-universal lat- 
tice scale physics (A ~ cAe -27r / 9 for small g). However, 
the long distance physics depends on these lattice scale 
effects only through the value of A, and has no direct 
dependence on g or A. 

We shall study correlators of the magnetization den- 
sity, M a (x,T) = SA/SH a (x,r). In the Hamiltonian for- 
malism, this magnetization is measured by the operator 
M a (x), and we shall focus on the real time, finite T cor- 
relation function 



C a p(x, t) 



AHt 1 



M a (x)e- tHt Mp(0))- (M a }(Mp 



where TL is the Hamiltonian corresponding to the action 
A, and the expectation values are with respect to the 
density matrix e~ n ^ T /Tre~ n ^ T . The dimensions of M 
are inverse length, and because M is a conserved den- 
sity, it does not acquire any anomalous dimension (i.e. 
no prefactors of powers of A or In A are required to ob- 
tain a finite A — > 00 limit), and its correlators are simply 
universal functions of combinations of x, t, c, T, H and 
A which are consistent with naive dimensional analysis 
in lengths and timestHl. For H <C A (which we assume 
throughout), these correlators describe the crossover be- 
tween two distinct limiting physical regimes: {%) T <C A, 
the 'quantum-disordered' regime, where strong quantum 
fluctuations create a paramagnetic ground state, and the 
excitations consist of a triplet of particles with energy 
(A 2 + c 2 ^ 2 ) 1 / 2 at momentum p, and («) A <C T <C cA, 
the high T regime of the continuum theory, where quan- 
tum fluctuations are marginally subdominantlij (by a fac- 
tor of l/ln(T/A)), and the excitations are a doublet of 
spin-waves about a locally ordered state; however ther- 
mal fluctuations of classically interacting spin-waves lead 
again to a paramagnetic state. The crossover between 
these regimes has been described for the static_-iiniform, 
spin susceptibility, Xi of the 0(N — 00) modelEJ. 
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FIG. 1. Two particle collision described by the S'-matrix 
(0). The momenta before and after the collision are the same, 
so the figure also represents the spacetime trajectories of the 
particles. 

In this paper, we shall obtain the space-time dependent 
C(x, t) of the 0(3) model in the low T region T < A. 
The ratio MJT is however allowed to be arbitrary. A re- 
cent papered computed C(0, t) by arguing that the triplet 
of particles could be considered free at low enough T. Ac- 
tually, such an approach is valid only for \t\ shorter than 
the mean collision time ~ e A / T /T (see below), and it is 
essential to include particle collisions at longer \t\ to ob- 
tain the crossover to diffusive behavior. Our semiclassical 
approach does this, and is valid for all \t\ ^> 1/T. 

There are two key observations that allow our exact 
computation for T <C A. The firstM is that as there is 
an excitation gap, the density of particles ~ e~ A / T , and 
their mean spacing is much larger than their thermal de- 
Broglie wavelength ~ c/(AT) -1 / 2 ; as a result the parti- 
cles can be treated semiclassically. In particular, taking 
the field H pointing along the 3 direction, the density of 
a particle with longitudinal spin m (m = —1, 0, 1) is 

„ - f ^ „-(A-mfl+c 2 P 2 /2A)/T _ . f^K -(A-mH)/T 

and therefore the total density p = p-\ + po +p% , and the 
magnetization (M a ) = (pi — p_i)6 a 3- The second obser- 
vations is that collisions between these particles are de- 
scribed by their known two-particle S'-matrixQ, and only 
a simple limit of this S'-matrix is needed in the low T 
limit. The r.m.s. velocity of a thermally excited particle 
vt = c(T/A) 1 / 2 , and hence its 'rapidity' ~ vt/c <S, 1. In 
this limit, the S'-matrix for the process in Fig [j] isQ 

^m\,ra'^ = (~^)^m 1 rn' 2 ^m 2 rn' 1 - (1) 

In other words, the excitations behave like impenetrable 
particles which preserve their spin in a collision. Energy 
and momentum conservation in d — 1 require that these 
particles simply exchange momenta across a collision 
(Fig |). The (-1) factor in (0) can be interpreted as the 
phase-shift of repulsive scattering between slowly moving 
bosons in d = 1. Indeed, it appears that the simple form 
of (Q) is due to the slow motion of the particles, and is 
not a special feature of relativistic continuum theory: we 
conjecture that (|l|) also holds for lattice Heisenberg spin 
chains in the limit of vanishing velocities. 

We now evaluate C(x,t) along the lines of a recent 
computation for the Ising modeti. We represent C(x, t) 

as a 'double time' path integral, with the e~ %nt factor 




1 ► X 

FIG. 2. A typical set of particle trajectories contributing 
to C(x,t). Each trajectory represents paths moving both for- 
ward and backward in time. The particle co-ordinates are 
Xk(t), with the labels k chosen so that Xk(f) < xi(t) for all t 
and k < £. Shown on the trajectories are the values of the 
particle spins mk which are independent of t in the low T 
limit. 

generating trajectories that move forward in time, and 
the e lW * producing trajectories that move backward in 
time. In the classical limit, stationary phase is achieved 
when the trajectories are time-reversed pairs of classi- 
cal paths (Fig ^) . Each trajectory has a spin label 
which obeys (|l|) at each collision; however as each col- 
lision contributes both to the forward and backward tra- 
jectories, the net numerical factor is simply +1. All 
of this impliesliil that the lines in Fig || are indepen- 
dently distributed uniformly in space, and with an in- 
verse slope determined by the velocity v which is dis- 
tributed according to the classical Boltzmann probabil- 
ity density V(v) oc e~ Av / 2c T . The spin, m, is as- 
signed randomly at some initial time with probability 
f m ee p m /p = e mH / T /(l+2 cosh(F/T)), but then evolves 
in time as discussed above (Fig ||). 

We label the particles consecutively from left to right 
by an integer k (see the caption of Fig ^) ; then their spins 
mfc are independent of t, and we denote their trajectories 
Xk(t). The longitudinal correlation C33 is given by the 
correlators of the classical observable 

M 3 (x,t) = J2 m kS(x-x k (t)) (2) 
k 

in the classical ensemble defined above. Now because 
the spin and spatial co-ordinates are independently dis- 
tributed, the correlators of nik and Xk factorize. The 
correlators of the rrik are easily evaluated: 

(m k mt) = Ax + A 2 Ske (3) 

where A, = (/j - f^ f and A 2 = h + JU - (/j - f^ f 
are simple, dimensionless, known functions of H/T only. 
Using (||) we have 

C 33 (x-x',t-t')=A 1 ({p{x,t)p(x' ,t')) - p 2 ) 

+ A 2 J2(^-x k (t))6(x' -z fe (i'))) (4) 
fe 
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where p(x,t) — J2k ~ x k(t)) is the spacetime depen- 
dent total density, all averages are now with respect to 
the classical ensemble, and (p(x, t)) — p. The two-point 
correlators of p(x, t) are also easy to evaluate: if the spin 
labels are neglected, the trajectories in Fig|^ are straight 
lines, and the density correlators are simply those of a 
classical ideal gas of point particles. The second corre- 
lator in (^), multiplying A 2 , is more difficult: it involves 
the self correlation a given particle k, which follows a 
complicated trajectory in the way we have labeled the 
particles (e.g. the trajectory of the —1 in Fig ||). For- 
tunately, precisely thiSpCorrelator was considered thtec 
decades ago by Jepsenfij and a little later by othersEj; 
they showed that, at sufficiently long times, each such 
particle executes free Brownian motion. Inserting their 
results into ([|), we obtained the final results presented 
below after some straightforward simplifications. 

An important property of the results is that they can 
written in a 'reduced' scaling formll3 determined by the 
classical dynamics. From the many independent param- 
eters c, A, T, and H, only a single length (L x ) and a sin- 
gle time (L t ) scale controls their spacetime dependence. 
These scales can be chosen to be 



F, 



1 

P 



1 



A 



p \2c 2 T 



1/2 



(5) 



Notice L x ~ ce A ' T / V AT is the mean spacing between 
the particles , and L t ~ e A / T /T is a typical time between 
particle collisions as vt — L x /L t V2. Our final result is 



C 33 (x,t) =p 2 



Ai.Fi 



A0F0 



L x Lt 



(6) 



where p 2 F\ is the connected density correlator of a clas- 
sical ideal gas in d = 1, 

F 1 {x,f)=e-* 2 / P /tsfr, (7) 
and p 2 F 2 is the correlator of a given labeled particleElIll, 



F 2 {x,t) 



(2Gi(u)Gi(-u)+fi(x,t)) 

x I Q (2iy/G 2 {u)G 2 (-u)) 

Gi(u)G 2 {-u) + G\{-u)G 2 (u) 



^G 2 {u)G 2 {-u) 



(2iy/G 2 (u)G 2 (- 



-(G 2 («)+G 2 (-«))t 



(8) 



with u = x/t, Gi(u) = erfc(u)/2, and G 2 (u) — 
e~ u /(2y/n) — uG\(u). These expressions satisfy 
J °° dxF-L <2 (x, t) = 1/2, which ensures the conserva- 
tion of the total magnetization density with time. For 
\t\ <C \x\ <C 1, the function F 2 has the ballistic form 
F 2 (x,i) w Fi(x,t), while for \t\ » 1, \x\ it crosses over to 
the diffusive form 



F 2 (x,t) 



(47Ti 2 ) 1 /4 



for large t. 



(9) 



In the original dimensionful units, (^) and (|9|) imply a 
spin diffusion constant, D s , given exactly by 



D„ = 



c 2 e A / T 



A(l + 2cosh(i?/T))' 



(10) 



Let us now consider correlations of the transverse mag- 
netization. It is convenient to work with the circularly 
polarized components of the magnetization M± = Mi ± 
iM 2 . The analog of @ is now M± — J^k S±kS(x—Xk(t)), 
where S + k {S-k) are the spin raising (lowering) opera- 
tors of particle k. In the double time path integral for 

C \-(x,t) the spin on some particle k is raised at time 

t = Q in the forward trajectory; at time t the lower- 
ing operator must act on the same particle or otherwise 
the trace over the classical trajectories vanishes. No- 
tice also that there is no raising or lowering of spins 
in the backward trajectory. As a result, the path in- 
tegral picks up a factor of e lHt from the additional Berry 
phase accumulated during the time the spin is raised dur- 
ing the forward trajectory, which is not compensated by 
the backward trajectory. This phase is multiplied by the 
self-correlation of the particle whose spin was raised, a 
quantity we have obtained above. Similar considerations 
apply to C-| and the final results are 



C T ±(x,t) = 2p 2 e ±lHt A T F 2 



Lt 



(11) 



where A T = f + f Tl . 

Next, we compute the local dynamic structure fac- 
tor S a p(uj) = f^dte-^Capfat). A subtlety arises m 
computing this Fourier transform. Notice that at short 
\t\ <C Lt, we have the ballistic behavior C(0, i) ~ 1/1*1) 
and so the t integral is logarithmically divergent. Our 
semiclassical results are valid only for \t\ ^> 1/T, and 
so we should cut-off the integral at small t, leading to a 
contribution ~ ln(bT/uj) where b is a numerical factor of 
order unity. In fact, it is possible to determine b precisely: 
at these-short times the earlier free quantum particle ap- 
proacbli3 is valid, and we determine b by matching the 
logarithm to their results. In physical terms, the short 
time cut-off is provided by the wavelike nature of the 
individual particles, at a scale where collisions are unim- 
portant. Our final results for S(uj) are 



S 33 M= - c J ^ [Ax {HTLt) + *i(V5F|«|i t )} 




-A 2 {ln(TL t ) + ^ 2 (V^HL t )}] 



{ln(TL t ) + <f> 2 (V^\uJTH\Lt)} 



The \n(TL t ) terms logarithmically violate the purely 
classical, reduced scaling formal^, and were fixed by 
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FIG. 3. The crossover function $2(fi) appearing in the lo- 
cal dynamic structure factor. 



matching to the short-time quantum calculation!! 4 ], 
scaling functions $i,2(^) were determined to be 



The 



)= In 



*i(n) 
$ 2 (n)=$i(n) 



7r[(V4 + fi 2 + 2) 1 / 2 - 



In 



4v^(V4 + fi 2 + 2) 1 /2 
(l + f2 2 /* 2 (r2)) 1 /2(i + ^(f7)) 
217 



(12) 



where 7 = 0.57721 ... is Euler's constant, and ^(fl) = 
{fly/1 + tt 2 /4 - fl 2 /2) 1 / 2 . We show a plot of the scal- 
ing function $2(0) in Fig ||: it clearly shows the ex- 
pected crossover from the large frequency ballistic be- 
havior $2(^ - ¥ 00) = ln(l/f2), to the small frequency 
diffusive form <J> 2 (ft -> 0) = 7r/(2 v / fi). 

The longitudinal relaxation rate of nuclei coupled to 

the electronic spins is 1/Ti = (Y/2)S^ (wjv), where T is 

determined by the electron-nucleus hyperfine coupling, 
and Wat is a nuclear frequency which can safely be set to 
zero. It is useful to explicitly note the H -C L^ 1 limit of 



1/Ti, where from (|12|), we have 



1 



TT X rAe~ 3A / 2T 



y/2D s H 



3T 



-A/T 




where x = ^ m H^o(pi ~ P—i)/H was known earlier 
For experimental comparisons, an important property of 
the above, pointed out to us by M. Takigawa, is that 
the low T activation gaps for 1/Ti (A 1 / Tl ) and x (A x ) 
satisfy Ai/ Tl /A x = 3/2. 

A quantitative comparison of our results with experi- 
ments requires a detailed study of the H and T depen- 
dencies of 1/Ti, along with consideration of effects due 
to spin-anisotropies and inter-chain couplings which can 
become important at low T and H. Such an analysis will 
be presented elsewhere; here, we simply note some trends 
which appear to receive a natural explanation from our 
theory. Values for the activation gaps A 1/7^ and A x 
been quoted for a number of experimental systemsEf 
and it has consistently been found that A 1 / Tl is larger 
than A x . In the spin 5 = 1 chain compound AgVP2Se 
Takigawa et. alB estimated A 1 / Tl /A x = 1.3; for the 



spin 5=1 chain compound Y^-BaiVzOs, Shimizu et. a/.i 
measured A\/t 1 /A x = 1.53 ± 0.08; finally, in the two- 
leg 5=1/2 ladder compound SrCuzOs, Azuma et. al.u 
found Ai/ Tl /A x = 1.6. 

Takigawa et. aZ.B also observed the diffusive 1 / \fH 
dependence of 1/Ti, from which the value of D s was es- 
timated: D s /a 2 fa 5.5 x 10 14 sec" 1 at T = 2%0K, where a 
is the lattice spacing. From measurementaill of x we may 
obtain A = 320-ftT , and c/a = 3.32A, which when inserted 
into (fil]) give D s /a 2 = 6.6 x 10 14 sec^, 1 . However, it 



ifii 



should be noted that numerical analysia_3 on the nearest- 
neighbor 5=1 antiferromagnet gives c/a — 6. 06 A, but 
using this value of c would also lead to a discrepancy in 
the theoretical prediction for x- 1—1 

Finally, we note that similar methodstil can be used to 
obtain dynamic, T > 0, correlators of the n a field: this 
will be described elsewhere. 

We are indebted to M. Takigawa for surveying and in- 
terpreting the experimental situation for us. We thank 
Satya N. Majumdar and T. Senthil for stimulating dis- 
cussions which provoked our interest in this problem, and 
I. Affleck for helpful remarks. This research was sup- 
ported by NSF Grant No DMR 96-23181. 



1 M. Takigawa et. al, Phys. Rev. Lett. 76, 2173 (1996). 

2 T. Shimizu et. al, Phys. Rev. B 52, R9835 (1995). 

3 M. Azuma et. al. Phys. Rev. Lett. 73, 3463 (1994). 

4 F.D.M. Haldane Phys. Lett. 93A, 464 (1983). 

5 E.S. Sorensen and I. Affleck, Phys. Rev. B 49, 13235 
(1994). 

6 A.M. Polyakov, Phys. Lett. B 59, 87 (1975). 

7 A.B. Zamalodchikov and A.B. Zamalodchikov, Ann. of 
Phys. 120, 253 (1979). 

8 Form factors in completely integrable models of quantum 
field theory by F.A. Smirnov, World Scientific, Singapore 

(1992) . 

9 A.M. Tsvelik, Zh. Eksp. Teor. Fiz. 93, 385 (1987) (Sov. 
Phys. JETP 66, 221 (1987)). 

10 J. Sagi and I. Affleck, Phys. Rev. B 5 3, 9188 (1996). 

11 S. Sachdev and A. P. Young, preprint cond-mat/9609185. 

12 A.V. Chubukov, S. Sachdev, and J. Ye, Phys. Rev. B 49, 
11919 (1994). 

13 M. Luscher, Phys. Lett. 118B, 391 (1982). 

14 Th. Joliceour and O. Golinelli, Phys. Rev. B 50, 9265 
(1994). 

15 D.W. Jepsen, J. Math. Phys. 6, 405 (1965). 

16 J.L. Lebowitz and J.K. Percus, Phys. Rev. 188, 487 (1969). 

17 M. Takigawa et. al, Phys. Rev. B 52, R13087 (1995). 

18 E.S. Sorensen and I. Affleck, Phys. Rev. Lett. 71, 1633 

(1993) . 



4 



